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The spectral branching behavior of the 2x2 operator matrix of the magneto-hydro- 
dynamic a 2 — dynamo is analyzed numerically. Some qualitative aspects of level crossings 
are briefly discussed with the help of a simple toy model which is based on a Z2— graded- 
pseudo-Hermitian 2x2 matrix. The considered issues comprise: the underlying SU(1, 1) 
symmetry and the Krein space structure of the system, exceptional points of branching 
type and diabolic points, as well as the algebraic and geometric multiplicity of correspond- 
ing degenerate eigenvalues. 
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1 The 2 x 2 operator matrix of the MHD a 2 dynamo 

The magnetic fields of planets, stars and galaxies are maintained by homoge- 
neous dynamo effects, which can be successfully described within magnetohydrody- 
namics (MHD) 1 . This is achieved by appropriately combining the Maxwell equa- 
tions of electrodynamics with the Navier-Stokes equations of hydrodynamics — for 
certain topologically non-trivial, helical flow fields. The resulting highly complicated 
equation systems are subject of large and cost-intensive computer simulations (see, 
e.g., Ref. Recently, the homogeneous dynamo effect has been demonstrated in 
large scale liquid sodium experiments in Riga/Latvia 3 and Karlsruhe/Germany 
pi]. Next generation experiments are currently planned at 7 sites around the world. 

One of the simplest dynamo toy models, which can be regarded as similar impor- 
tant for MHD dynamo theory like the harmonic oscillator for quantum mechanics, 
is the spherically symmetric a 2 — dynamo pQ in its kinematic regime. Its operator 
matrix has the form [5] 

-m) w 

and consists of formally selfadjoint blocks 

Q[a] :=pap + a l -^±}±. (2) 

*) Presented at the 2nd International Workshop "Pseudo-Hermitian Hamiltonians in Quantum 
Physics", Prague, Czech Republic, June 14-16, 2004 
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(p = —i(d r + 1/r) is the radial momentum operator.) It describes the coupled 
I— modes of the poloidal and toroidal magnetic field components in a mean-field 
dynamo model with helical turbulence function (a— profile) a(r). The differential 
expression Q has the symmetry property [5] 

= jtf/MJ, J=(J J) (3) 

and shows that — depending on its domain T>(Hi[ce\) and the chosen boundary 
conditions for the two-component eigenfunctions ifi — the corresponding operator 
can be pseudo-Hermitian (J— self adjoint). Modulo J— selfadjoint extensions, this 
will hold for functions ip with idealized boundary conditions at r = 1 

i . „/, a — n-/ en n-/ n-/ — t . fn 



V(Hi[o\) := j^> = ^ J : ^ ^n = H®H, H = L 2 {fl, r 2 dr), 

n=[0,l], V(1) = 0, n/>MU o -+0}, (4) 
whereas the physically realistic boundary conditions 

-B;V =0, Bi= diag[0 r + + l)/r, 1] (5) 

r— 1 

lead to a non— J— Hermitian operator. 

The dynamo effect starts when the first eigenvalue A of Hi [a] enters the right 
half-plane 5ft(A) > 0. This leads to an exponential growth ~ e At of the solutions of 
the time dependent kinematic dynamo problem — and a corresponding growth of 
the magnetic field components. In Ref. |S] it was demonstrated numerically that 
this can happen even on a complex-valued branch of the spectrum, so that the 
dynamo can start in an oscillating regime. 

In figures n - El some typical spectral branches are shown for an a 2 — dynamo 
operator with physically realistic boundary conditions (j^J . The real and imaginary 
components of the eigenvalues for modes with angular mode number I = 1 and 
radial mode numbers n = 1, ... ,9 are depicted over a scaling parameter C of an 
a(r)— profile which is chosen as quartic polynomial a(r) = C x [ag + a,2r 2 + a^r 3 + 
0,4 r 4 ]. Of special physical interest are the critical value C c , where the spectrum 
enters the right half-plane 3?(A) > 0, as well as the location of those level crossing 
points, where two real-valued branches (non-oscillating regime) of the spectrum 
meet and continue to evolve as two complex conjugate branches (oscillating dy- 
namo regime). Within the considered range of C such transitions occur locally only 
pairwise — there are always only two spectral branches which are locally involved 
in such transitions. As shown in fig. 01 globally there are more branches involved 
in mutual transitions. A detailed and rigorous non-numerical study is still missing. 
Below, we collect some few building blocks for such an analysis. 

Some first and rough qualitative aspects of the behavior of the dynamo system at 
a level crossing point can be easily understood, e.g., by passing from the eigenvalue 
problem for the linear pencil 

Li[a,\]i/>:=(6t[a]-\)il) = (6) 
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Fig. 1. Large-scale behavior of real and imaginary parts K(A), Sr(A) of the a 2 — dynamo 
spectrum (for physically realistic boundary conditions |0) as function of the scaling pa- 
rameter C in the concrete a— profile 

a(r) = C x (l - 26.09 x r 2 + 53.64 x r 3 - 28.22 x r 4 ). 



of the 2 x 2— operator matrix via substitution 

*-(iB(D + *!*)■ ° w ^° (7) 
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to the equivalent eigenvalue problem of the associated quadratic operator pencil 

L t [a,X]iPi = {[Q[l] + A]^[Q[1] + A]-Q[a]}v>i=0 

= (A 2 \ 2 + A 1 \ + A )yj 1 = 0. (8) 
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Fig. 3. Zooming into another scaling region with multiple level crossings. Interestingly, 
although locally only two branches undergo a transition from pairwise real- valued eigen- 
values to pairwise complex-conjugate ones, globally such transitions occur also for more 
branches. In the depicted scaling region three branches participate in mutual transitions. 
This is a natural indication of the underlying Riemann surface structure of the operator 

spectrum (see, e.g., 003)- 



Solving the quadratic (functional) equation 

Mi[a, A] := (L t [a, X]ipi, tpi) = a 2 A 2 + aiX + a = 0, aj := (Ajipi,ipi) (9) 
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Fig. 4. The real and imaginary parts of crossing spectral branches without transition from 
real to complex eigenvalues "feel" each other. 



for A we see from the solutions 



A ± = — y-ai ± \ja{ - 4a a 2 J 



(10) 



that a level crossing occurs when the discriminant A = a\ — 4ao<22 of Eq. 
vanishes. At the level crossing points it holds 



M l [a,X]=0, d x M l [a,X}=0 



(11) 
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Fig. 5. The long-term parallel location (repelling) of different branches without real-to- 
complex transition is not a numerical artefact. The branches can be faithfully identified 

as different ones. 



and the operator pencil has a Jordan-Keldysh chain |§1 111)1 ITT] {ipi, Xi> <^i}j con ~ 
sisting of the eigenfunction ipi and the associated functions xi, 4>i which satisfy the 
relations 

Ll[a,Ao]^ = 0, (12) 
L l [a,X \xi + d x Li[a,X]\ Xo il}i = 0, (13) 
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Li[a,X )(j) 1 + d x Li[a 1 \]\ Xa xi + -dlLi[a 1 \]\ XQ ^ 1 = 0. (14) 

Due to the lack of rigorous analytical and non-numerical operator theoretic studies 
of the oc — dynamo, we will collect in the next section some few qualitative facts 
about a highly simplified 2x2 matrix model with some rough structural analogies 
to the operator matrix (JTJ. 

2 Z2 graded pseudo-Hermiticity 

In this section 1 ) we start from the simple model of an a 2 — dynamo with idealized 
boundary conditions Q so that the corresponding operator Hi[a] is J— selfadjoint 
( J— pscudo-Hermitian) [see relation ©]. We will try to get a rough intuitive insight 
into some of the basic properties of such a system in the vicinity of a level crossing 
point. 2 ) 

General pseudo-Hermitian operators H were defined in Refs. as operators 
which satisfy a relation 

H = ijH^r 1 - 1 . (15) 

Because of the involution property of J, i.e. J 2 = I, J = J -1 , the operator Hi[a] 
is part of a narrow subclass of operators with H = \xH* ^l, /i 2 = /, /i -1 = /i. To the 
same subclass belong the VT— symmetric Hamiltonians of Refs. |151 1161 IT71 118j . 
which are V— pseudo-Hermitian and for which necessarily holds V 2 — I, as well as 
the Schrodinger Hamiltonian associated to a Wheeler-DeWitt equation of minisu- 
perspace quantum cosmology considered in Ref. |14j . A general feature of systems 
with involutive pseudo-Hermiticity operators [i is an underlying Z2— graded struc- 
ture of their Hilbert space: Every state (vector) of the Hilbcrt space can be naturally 
split into (projected onto) /x— even and \i— odd components (for VT— symmetric 
Hamiltonians into parity-even and parity-odd states, see, e.g., |T3IIT5] ) 

x = P+x + P-x = x + + x_, P±:=-(I±(i) (16) 

with 

fJLX± = ±x±. (17) 

Explicitly, this leads to a natural Z2— grading of the Hilbcrt space H H+ © Tt- 
and a Krein space structure 3 ) |21|: Beside the usual inner product (., .) of the Hilbert 
space TL, which induces a non-negative norm (x,x) = {x+,x + ) + (a;_,a;_) > for 
x± 6 H±, one can consider a Krein space [., .] M ) with indefinite inner product 

1 ) The few issues presented in this and the next section are part of a detailed study given in: 
U. Giinther, Z2— graded pseudo-Hermitian systems: exceptional points, pseudo-unitary fibrations 
and nontrivial holonomy, in preparation. 

2 ) After finishing this proceedings contribution we became aware that some aspects of level 
crossings in Z2— graded pseudo-Hermitian systems had been briefly discussed earlier — with the 
help of different techniques — in Refs. |12J|13J. The present analysis overlaps only marginally with 
those results and provides a different view on the subject. 

3 ) For a detailed introduction into the operator theory over Krein spaces see, e.g., Refs. II 1II2UI . 
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[x, y]^ = {nx,y) = (x+,y+) — (x— (and a corresponding indefinite "norm"). 
One naturally distinguishes states (vectors) of positive type [a;,a;] M > 0, of negative 
type [x, x] M < 0, and isotropic states [x,x] M = 0. (In rough analogy to Minkowski 
space this corresponds to time-like, space-like, and light-like vectors.) Because of 
fiH = W /i a fi— selfadjoint (ji— pseudo-Hermitian) operator is selfadjoint in the 
Krein space 

[Hx, y}^ = {fiHx, y) = (x, fiy) = (x, fiHy) = [x, Hy}^ . (18) 

A natural representation of a Z2— graded system can be given in terms of 2- 
component vectors and 2 x 2— operator matrices 




where H ++ = P + HP + , etc. The fi— pseudo-Hermiticity implies 

H++ = Hl + , ff__ = ffl_, ff+_ = -ffi + . (20) 

The operator if; [a] of the a 2 — dynamo with idealized boundary conditions can 
be transformed into this representation by first diagonalizing the involution opera- 
tor J 

J^V = S- 1 JS, S=^(\ (21) 

Applying then the same transformation S to Hi [a] and the elements of the Hilbert 
space H one obtains the equivalent operator [2] 

H l [a]=S-"H l {a]S= 1 -( Q ^- 2]+a n 7 Q[a] ^ ) (22) 
L J L J 2 V Q[a] - a Q[-a - 2] - a J K ' 

which acts on 2-vectors 



( 0+ " 




lp2 


+ ^1 \ 






v "02 


-^1 J 



(23) 



We leave a detailed analysis of this highly non-trivial operator to future studies. 

Instead we try to get a rough qualitative understanding of the level-crossing in 
a Z2— graded system. For this purpose, we consider its simplest example 4 ) — the 

4 ) This gives a level crossing with transition from a pair of real- valued eigenvalues to a pair of 
complex conjugate eigenvalues. Apart from such level crossing with real-to-complex transitions, 
there occur level crossings without such transitions. They are present in all figs. 13 and can 
be easily understood in a matrix setup as crossings of spectral branches which belong to different 
Z2— blocks. 
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eigenvalue crossing of a 2 x 2 matrix with the basic symmetry properties 120(1 : 

«=(-*• !)-(? l) +k ' -f 

e = (a + d)/2, f = (a-d)/2, b = 6i + %, a, d, &i, 62 € R. (24) 
For the eigenvalues holds 



det [H - EI) = £ = e ± - 6? - b\ (25) 

and we see that a level crossing occurs at 

A(f,b 1 ,b 2 ):=f 2 -b 2 1 -b 2 2 = 0. (26) 

A(/, 61,62) = defines a two-dimensional variety (a double-cone) in the three 
dimensional parameter space M. 3 (/, 61, 62) so that for / ^ ^ |6| the degeneracy 
has co-dimension one. This is in obvious contrast to a Hermitian 2x2 (spin) matrix 
problem with 

hs = ( I b f ) (27) 

and eigenvalues 

det {H s - EJ) = £ s = e ±^/ 2 + 6 2 + 6 2 . (28) 

Here a degeneracy A s (/, 61, 6 2 ) := f 2 +b\+b\ = occurs only in the single (diabolic) 
point / = bi = b 2 = and the crossing has co-dimension three [221 1221 ■ 

Another difference between the two models is the underlying symmetry of "iso- 
energetic" (adiabatic) deformations. Obviously, the eigenvalues E and E s in (|25|l 
and (|28fl are invariant, respectively, under SO(l,2) and SO(3) transformations in 
the parameters /, 61 , 6 2 . This is in natural correspondence with the generators 
of " iso-energetic" (adiabatic) transformations which are defined by the matrices 
H and H s themselves. Modulo the Abelean U(l) transformations induced by the 

elements ieo ( jt !?) these generators are Lie algebra elements of the type 



v 1 , 

ih = -&20-1 - &i02 + G su(l, 1) ~ so (l,2) ~si(2,R), (29) 
«ft. s = i (pi&i + 6 2 (j 2 + 63C3) g su(2) <~ so(3). (30) 

((Tj are the Pauli matrices.) Of course, this interlinks 5 ) also nicely with the invari- 
ance transformations of the two-dimensional Krein space 6 ) (J^ , [., .] M ), \i — (T3 and 

5 ) Because of its underlying pseudo- unitary symmetry pseudo-Hermitian ("PT— symmetric) 
quantum mechanics is not a special variant of quaternionic quantum mechanics 1251 . One sees this 
immediately from the simplest models with H and H a as Hamiltonians. Quaternionic structures 
are naturally connected with SU(2) symmetries 1241 1251 and H a , whereas the pseudo-Hermitian 
model with Hamiltonian H has invariants which are related to SU(l, 1) symmetry transforma- 
tions. Additionally, we note that SU(2) is a compact group whereas SU(1, 1) is a non-compact 
one. 

6 ) Because of dim(T-i-i-) = 1 < 00 for the considered 2x2 matrix model, this Krein space is a 
so called Pontryagin space II 1| . 
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the 2-vector Hilbert space 7i s of the spin system. The inner product [., is in- 
variant under pseudo-unitary {7(1,1) ~ {7(1) x S{7(1, 1) transformations, whereas 
the inner product (., .) in 7i s is invariant under unitary {7(2) ~ {7(1) x SU(2) 
transformations. 

Further insight into the structure of the simple pseudo-Hermitian 2x2 matrix 
eigenvalue problem can be gained by diagonalizing the matrix H. Before we do this 
explicitly, we note that any 2x2 involution matrix 77 

rf = I, 77 = r^ 1 = (31) 

belongs to one of the following classes 77+ or 7i_ 

det(77+) = l, 77+ e {/,-/}, (32) 

det(r?_) = -1, r\- = a\cr\ + a 2 (T 2 + a 3 <7 3 , a 1 +a 2 + a 3 = l. (33) 

This means that two involution matrices Tji-, 772- are connected by an ST/ (2) ~ 
Z 2 x 50(3) rotation. 

For non-degenerate eigenvalues E the diagonalization is done as 



H = SDS~\ D= 60 " 6 u , . 5= ~ 71 - 
V e + e / V 7l 72 

e := A 1 / 2 , 7li2 G C. (34) 



The complex constants 71,2 are still arbitrary and we can fix them by requiring 
that the diagonalization preserves the pseudo-Hermitian structure, i.e. that there 
exists an involution matrix 77 so that 

D = rjD^rj. (35) 



From the substitution chain 

H = nH*ii =>• S77S" 1 =M(S" 1 ) t i7 t SV D = S'V (S" 1 ) 1 D^S^S (36) 
and relation l|35l) one identifies 

rj = S^S = S- 1 ^i(S- i y = fr 1 = ^ (37) 

so that 

M = SS t /j5S t , rj = S^Sr]S^S. (38) 

From the latter equations one concludes that | det(S)| 2 = 1. Explicit calculations 7 ) 
show that one can set det(S) = 1 and that S is in general neither unitary, 5* 7^ S _1 , 
nor pseudo-unitary. Complementary information can be gained from the explicit 



7 ) See footnote 111 
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structure of the matrix D. It holds 

A>0:=^ e = A 1 /2>o, D = ( 6 ° ~ e ° ) = (39) 



A<0:=^ e = «|e|, D 



e — i|e| 
e + i|e| 



e + i|e| 
e -i|e| 

and one obtains from l|35l) and the explicit form of S in H34[) (by appropriately 
tuning the constants 71,2) that 

A > : / > : ry = ( J \ , (41) 

/<0: V = ( J _°j ), (42) 

A<0: „=(j J). (43) 

We see that smooth changes of the starting fi— pseudo-Hcrmitian matrix H lead to 
qualitative " switchings" (discontinuities in the mapping /x 1— » 77 which correspond to 
rotations in the space of 2 x 2 involution matrices r\- — see Eq. H330 1 in the pseudo- 
Hermitian structure of the diagonal matrix D. These "switchings" occur when the 
system intersects critical surfaces in the parameter space: (1) the surface / = 
and (2) the degeneration double cone A = 0. In the latter case, the "switching" 
in the pseudo-Hermiticity matrices corresponds to a " switching" from a Hermitian 
diagonal matrix D = for A > to a complex-valued diagonal matrix D 7^ 
for A < 0. A "switching" occurs also in the properties of the eigenvectors of the 
diagonal matrix D. Choosing these eigenvectors in the simplest form as 

D|±>=(e ±e)|±>, !+>=(?). I">=(J 

<±|±>=1, <+| = (0,l), < -I = (1,0), (44) 

one finds for the Krein space inner products 

A>0: [±,±], = ±Bign(/) ) [±,=F], = (45) 

A<0: [±,±] v = 0, [±,T] V = 1. (46) 

This means that the eigenvectors |± > are of positive or negative type for real- 
valued eigenvalues, A > 0, and of isotropic type for pair-wise complex conjugate 
eigenvalues, A < 0. We explicitly reproduced a basic result of Krein space theory 
which is discussed in |2U] and which is also implicitly present, e.g., in Refs. |14II17I 

us]. 
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3 Exceptional points 

From the structure of S in l|34l) one sees that in the degeneration limit e — > 
the determinant det(S) vanishes, det(S) — > 0, S 1 becomes singular and the diago- 
nalization i|34[l breaks down. Instead the matrix £) turns into a Jordan block. For 
e = 0, f 2 — \b\ 2 7^ 0, / = ±|6|, a = d± 2|6| one finds the explicit relation 

H = SDS- 1 , D =(q e)> S= ( T F t)' E = d±\b\. (47) 

Additionally, one observes 8 ) that only |— > survives as an (geometric) eigenvector 
of D, whereas |+ > is now an associated vector (algebraic eigenvector) 

(D-EI)\->=0, (D-EI)\+ >={->, (D- Elf\+ >= 0. (48) 

Hence, the degenerate eigenvalues on the double cone have algebraic multiplicity 
two and geometric multiplicity one. The degeneracy is a branching point degeneracy 
|77] — a double cone of exceptional points of branching type in the sense of 
Kato 

A singularity of higher order occurs at the center of the double cone, i.e. in the 
diabolic point [301 at the origin / = \b\ = of the parameter space M 3 (/, b 1: b 2 ). 
There the transformation matrix S in (|47|) becomes singular and instead of the 
Jordan block (|47|l the matrix H has the form 

H = ( f E = a = d. (49) 

This is the same type of co-dimension three degeneracy as in the case of the spin 
matrix H s — with two (geometric) eigenvectors, i.e. the algebraic and the geometric 
multiplicity of this eigenvalue coincide and equal two. 

The question of whether nontrivial holonomy (geometric phases |23|) and chiral- 
ity properties of cigcnfunctions in the vicinity of level-crossing points can play 
a role in the dynamics of unstable magnetic field configurations with field reversals 
|2] remains as one of the interesting open issues. 

We thank the Professors M.V. Berry, W.D. Heiss, H. Langer, C. Tretter and M. Znojil 
for useful comments at different stages of the work. U.G. acknowledges support from DFG 
grant GE 682/12-2. 



8 ) The existence of a Jordan block structure with a single geometric eigenvector is a generic fea- 
ture of systems at level-crossing points and was observed for a ID pseudo-Hermitian Hamiltonian, 
e.g., in 26. 
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